g-Breathers and the Fermi-Pasta-Ulam Problem 
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The Fermi-Pasta-Ulam (FPU) paradox consists of the nonequipartition of energy among normal 
modes of a weakly anharmonic atomic chain model. In the harmonic limit each normal mode 
corresponds to a periodic orbit in phase space and is characterized by its wave number q. We continue 
normal modes from the harmonic limit into the FPU parameter regime and obtain persistence of 
these periodic orbits, termed here g-Breathers (QB). They are characterized by time periodicity, 
exponential localization in the g-space of normal modes and linear stability up to a size-dependent 
threshold amplitude. Trajectories computed in the original FPU setting are perturbations around 
these exact QB solutions. The QB concept is applicable to other nonlinear lattices as well. 

PACS numbers: 63.20.Pw, 63.20.Ry, 05.45.-a 



Fifty years ago E. Fermi, J. Pasta, and S. Ulam (FPU) 
published their celebrated report on thermalization of ar- 
rays of particles connected by weakly nonlinear strings 
bringing forth a fundamental physical and mathe- 
matical problem of the energy equipartition and ergodic- 
ity in nonlinear systems. Series of numerical simulations 
showed that energy, initially placed in a low frequency 
normal mode of the linear problem with a frequency uj q 
and a corresponding wave number q, stayed almost com- 
pletely locked within a few neighbour modes, instead of 
being distributed among all modes of the system. More- 
over, recurrence of energy to the originally excited mode 
was observed. 

Much effort has since been expended to understand 
and explain the FPU results (see 0- 0- 01 for reviews). 
Two major approaches were developed. The first one, 
taken by N. Zabusky and M. Kruskal, was to analyze 
dynamics of the nonlinear string in the continuum limit, 
which led to a pioneering observation of solitary waves 
0. The second approach, followed by F. Izrailev and 
B. Chirikov pointed to the existence of a "stochasticity 
threshold" in the original FPU system Q. For strong 
nonlinearity (or simply large energies) the overlap of non- 
linear resonances 0] leads to strong dynamical chaos, de- 
stroying the FPU recurrence and insuring fast conver- 
gence to thermal equilibrium. If nonlinearity is below a 
(size-dependent) "stochasticity threshold", the dynamics 
of the chain remains similar to that of the unperturbed 
(linear) system for large time scales. Later studies 0, 0| 
showed that the "local" dynamics of four consecutive low- 
frequency modes may become substantially chaotic, while 
almost all initial energy stays localized in these modes 
during the time of computation. The redistributed mode 
energies fall exponentially with increasing mode numbers 
in this regime (coined "weak chaos") and the energy flow 
to higher frequency modes was argued to be exponen- 
tially slow (due to Arnold diffusion). 

The results obtained through this second approach 



lead one to formulate some important questions. Firstly, 
since the dynamics below the stochasticity threshold is 
localized in g-space for long times, do time-periodic tra- 
jectories with almost all energy locked in a single mode 
for infinite times, coined q-Breathers (QB), exist, and are 
they close to the ones studied by FPU? Secondly, are the 
stability thresholds of such QBs related to the various 
stochasticity thresholds mentioned above? And finally, 
is the concept of QBs applicable also for other spatially 
extended nonlinear lattices, including generalizations to 
higher lattice dimensions? A strong motivation for this 
study is the fact that in the q-space representation we 
deal with oscillators which are uncoupled in the limit of 
small amplitudes, and, moreover, with frequencies being 
different for each mode. Nonlinearity induces coupling 
between oscillators. That is reminiscent of the case of dis- 
crete breathers (DB), that are time-periodic and spatially 
localized excitations on networks of interacting identical 
anharmonic oscillators, which survive continuation from 
the trivial limit of zero coupling Notably, DBs exist 
also in FPU lattices [ill ]. 

The FPU system is a chain of N equal masses cou- 
pled by nonlinear strings with the equations of motion 
containing quadratic (the a-model) 
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or cubic (the /3-model) 
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interaction terms, where x n is the displacement of the n- 
th particle from its original position, and fixed boundary 
conditions are taken xq = xn + i = 0. A canonical trans- 
formation x n (t) — J X] Qq(t) sin ( j^pj ) takes into 
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the reciprocal wavenumber space with N normal mode 
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coordinates Q q (t). The equations of motion then read 
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Qq + u\Q q 



for the FPU-a chain (QJ and 
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for the FPU- ,3 chain (0), where uj q = 2 sin (nq/2(N +1)) 
are the normal mode frequencies, and A q ^,j and C q .i.j iTn 
are coupling coefficients 0. For small amplitude ex- 
citations the nonlinear terms in the equations of mo- 
tion can be neglected, and according to © and <@J the 
(/-oscillators get decoupled, each conserving the energy 

E q = \ (Q q + ^qQq^j in time. Especially, we may con- 
sider the excitation of only one of the g-oscillators, i.e. 
E q ^ for q = qo only. Such excitations are trivial time- 
periodic and g-localized solutions (QBs) for a — [3 = 0. 

Let us consider the /?-model and choose first (3 = 0, 
excite a normal mode with q = qo to the energy E qa = E 
and let all other g-oscillators be at rest. With that we 
arrive at a unique periodic orbit in the phase space of the 
FPU model. In order to guarantee continuation of this 
periodic orbit into the case (3 ^ and following 0] we 
conclude, that it is enough to request the nonresonance 
condition nu) qo ^ ^ q ^ qo which is the generic case for a 
finite system size N (here n is an integer) (see also (l2^ . 
We expect that the orbit will stay localized in g-space at 
least up to some critical nonzero value of (3. The same 
argumentation can be applied to the a- model 

We compute QBs as well as their Floquet spectrum nu- 
merically using well developed computational tools [111 ], 
and compare the results with analytical predictions, de- 
rived by means of asymptotic expansions. As a zero-order 
approximation for the numerical computation we take the 

g -th linear mode: x n {t) = y jvq^Qg (*) sm (t^+t)- For 
the /3-model the initial conditions are Q q (t = 0) = 0, 
and Q qo (t = 0) = yj2E-Z& qo Q 2 q (t = 0). We map 

the space of y = {Q q } onto itself by integrating the ini- 
tial condition up to the time when Q qo {t) = again: 
y n+1 = J 7 (y n ). A periodic orbit is a fixed point of that 
map. The vector function Q = T(y) — y is used to calcu- 
late the Newton matrix Af = dQ(y)i/dyj. The iteration 
procedure y ' = y—J\f~ 1 G{y) continues until the required 
accuracy e is obtained: ||.F(27) — j7||/||j/1| < £ (we have 
varied e from 10 -5 to 10~ 8 ), where ||y|| = max{|j/i|}. 
For the a-model we used a modified scheme choosing 
x s (t = 0) = where s = [(N + l)/2go] corresponds to 
the antinode of the mode Q qo . We map the phase space 
r (excluding x s ) onto itself integrating until x s (t) = 
again. With the above notations we use a Gauss method 
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FIG. 1: Energy distributions between g-modes in QBs for dif- 
ferent nonlinear coupling coefficients (3 versus q in linear and 
log scales with analytical estimations of the QBs exponential 
localization (dashed lines). Parameters are E = 1.58, go = 3, 
N = 32. Only odd modes are shown (see text). The symbols 
for q 7^ 3, 9, 15, 21, 27 represent upper bounds, the real mode 
energies might be even less. Note that QBs persist even far 
beyond the stability threshold (see FigEJ. 
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FIG. 2: Absolute values of Floquet multipliers of QBs 
with the energy E = 1.58 and go = 3 and different N versus 
13. Symbols: numerical results, lines: analytical results. 



to solve the equations Q(r) = N(r — r ') for the new it- 
eration r ' and do final corrections to adjust the correct 
energy E. To compute the linear stability of the found 
QB, we linearize the phase space flow around it, and map 
that flow onto itself by integrating over one period of the 
QB. The corresponding symplectic Floquet matrix can be 
computed numerically and subsequently diagonalized. If 
all eigenvalues fi have absolute value one, the QB is sta- 
ble, otherwise it is unstable (ll] |. 

First, we apply our method to the /3-model with qo = 3 
and E = 1.58 which is very close to the value 1.5 chosen 
in We obtain QBs which are exponentially localized 
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FIG. 3: Stable QB solutions for a = 0.025 and a = 0.25, E = 
0.077, iV = 32, and go = 1; the latter corresponds to the 
original numerical FPU-a study Q. An analytical estimation 
of the QB exponential localization in case a = 0.025 is shown 
with a dashed line. 



in g-space (FigQ}. The smaller j3, the faster is the decay 
of the energy distribution with increasing wave number 
q. Note, that due to the parity symmetry of the /3-model 
(Eq.lQJ is invariant under x n — > —x n for all n) only odd 
q-modes are excited by the go = 3 mode and get coupled 
|l4j . In FigI2we plot the absolute values of the Floquet 
eigenvalues of the computed QBs versus for different 
system sizes N. QBs are stable for sufficiently weak non- 
linearities (all eigenvalues have absolute value 1). When 
j3 exceeds a certain threshold two eigenvalues get absolute 
values larger than unity (and, correspondingly, another 
two get absolute values less than unity) and a QB be- 
comes unstable. Remarkably, unstable QBs can be traced 
far beyond the stability threshold, and moreover, they re- 
tain their exponential localization in g-space (FigQJ. As 
(3 is increased further, new bifurcations of the same type 
are observed. 

By an asymptotic expansion of the solution to <@J in 
powers of the small parameter a = /3/(N+l) we estimate 
the shape of a QB solution Qi(t) localized in the mode 
go- The energies of the modes qo, 3go,. . . ,(2n + l)qo,. . . 
are given by 
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S^q 2 



(5) 



up to a an error of the order (2n + l) 2 q 2 /N 2 . Dashed 
lines in FigQ]are obtained using © and show very good 
agreement with the numerical results. 

Using standard secular perturbation techniques we ap- 
proximate the frequency f2 of the QB solution as ft = 

^qoV 1 ' RfJV+11 + w \( 



-)) . The instability threshold 



8(JV+1) T (N+l)' 2 > 

observed in Fig|2lcan be obtained analytically by making 
a replacement Qi = Qi (t) + £j in the equations of motion 
<@J and linearizing the resulting equations with respect 



£ + A£ + h(l + cos 2flt)B$, + 0(/i 2 )| = (6) 



is a vector, A = (SijUj 2 ) is a diago- 
= (bij) is a coupling matrix, and h — 



where £ = (&) 
nal matrix, B = 
3f3E/2(N + 1) is a small parameter. We analyze para- 
metric resonance in (jSJ , treating h and fl as independent 
parameters. In the limit h — > the equilibrium point 
£ = is stable for all values of O except for those which 
satisfy ujk + wj = 2nfl = 2nVl n ki where n > 1, and the 
modes k and I belong to the same connected component 
of the coupling graph whose connectivity is defined by 
the matrix B. 

We seek for a solution to © at £1 = Q, n ki{l +5), 6 = 
0(h), in the form £ = J2t=-oo f m e ( ~ lGl+z+2lm ^ t + c.c. 
where Co = 0^(1 + 5) — —u>i(l + S) + 2nfl, f m are un- 
known complex vector amplitudes, and z = 0(h) is a 
small unknown complex number. 

The nearest primary resonance corresponds to k = 
qo~ l,/ = go + l) n = 10|- In the vicinity of the bifurca- 
tion point the absolute values of the Floquet multipliers 
involved in the resonance are obtained as 



I A*ii j 



1 ± 



A(N + iy 



IR-1 + 



(7) 



where R = &f3E(N + 1)/tt 2 . The bifurcation occurs at 
R = l + 0(\/N 2 ). The result is plotted in FigElwith 
solid lines for N = 32, 64 and 128, demonstrating good 
agreement with the numerical results. The agreement 
improves with increasing The instability thresh- 

old for QB orbits (Fig|2j| which is obtained analytically 
using the parameter R iJJJ, coincides with the criterion 
of transition to weak chaos reported by De Luca et al @|. 

We have used one of the original parameter sets of the 
FPU-a study a = 0.25, E = 0.077, A - 32 Q (and add 
to that the case a = 0.025 as well for comparison) to find 
stable exponentially localized QB modes with most of the 
energy concentrated in go = 1 (FigEJ)- We use again an 
asymptotic expansion of the solution to (pj in powers of 
the small parameter p = a/^/2(N + 1) and obtain that 
the energies of the modes go, 2g ,. . . ,nq , are given by 



En 
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n A E„, 



a^Ei\N + lf/ 2 
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The dashed line in FigEI is obtained using © in case 
a = 0.025, E = 0.077, N = 32,q = I and shows very 
good agreement with the numerical results flf| . 

How are QB-solutions related to the original FPU stud- 
ies? A QB requires specific initial conditions, which were 
not used in earlier numerical studies. However, nearby 
quasiperiodic solutions are expected to retain stability 
and exponential localization of their QB generator for 
times long compared to the QB period. In Fig^lwe com- 
pare snapshots of displacements at different times ob- 
tained for the original FPU trajectory in |l| and for the 
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lattices and other Hamiltonians are straightforward, due 
to the weak constraint imposed by the nonresonance con- 
dition needed for continuation. QBs can be also expected 
to contribute to peculiar dynamical features of nonlinear 
lattices in thermal equilibrium, e.g. the anomalous heat 
conductivity in FPU lattices 

We thank T. Bountis, F. Izrailev, V. Shalfeev and V. 
Zakrzewski for stimulating discussions. M.I. and O.K. 
appreciate the warm hospitality of the Max Planck In- 
stitute for the Physics of Complex Systems. M.I. also 
acknowledges support of the "Dynasty" foundation. 



FIG. 4: Snapshots of displacements (a) of the original FPU 
trajectory for a = 0.25, E = 0.077, N = 32 Q and (b) of the 
corresponding exact QB solution from Fig|^|taken at different 
times. 

numerically exact QB solution from Fig for a = 0.25 
and observe similar evolution patterns. Moreover, we 
took a series of points on a line which connected initial 
conditions of the FPU trajectory (E q7 n = 0) with the 
numerically exact QB solution from FigG3 For each of 
these points we integrated the corresponding trajectory 
and measured the average deviation A from the QB orbit. 
The dependence of A on the line parameter turns out to 
be an almost linear one, starting from zero when being 
very close to the QB orbit, and ending with a maximum 
value when being close to the FPU trajectory. That sup- 
ports the expectation that the FPU trajectory is a per- 
turbation of the QB orbit. The FPU recurrence is grad- 
ually appearing with increasing A and is thus directly 
related to the regular motion of a slightly perturbed QB 
periodic orbit, which we tested also numerically. 

In conclusion, we report on the existence of g-breathers 
as exact time-periodic low-frequency solutions in the non- 
linear FPU system. These solutions are exponentially lo- 
calized in the q space of the normal modes and preserve 
stability for small enough nonlinearity. They continue 
from their trivial counterparts for zero nonlinearity at fi- 
nite energy (in that limit they simply correspond to one 
mode with wave number go being excited, and all the 
other modes being at rest). The stability threshold of 
QB solutions coincides with the weak chaos threshold in 
Persistence of exact stable QB modes is shown to be 
related to the FPU paradox. The FPU trajectories com- 
puted 50 years ago are perturbations of the exact QB 
orbits. Remarkably, localization in g-space persists even 
for parameters when the QBs turn unstable. The concept 
of stable QBs and their impact on the evolution of excita- 
tions in the FPU system is expected to apply far beyond 
the stability threshold of the QB solutions reported in 
the present work. Generalizations to higher dimensional 
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